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Abstract 

A  high-order  (10th)  mathematical  model  for  the  quiescent  geomag¬ 
netic  field  from  Epoch  1955.0  to  Epoch  1995.0  was  used  to  calculate 
the  asymptotic  directions  of  vertically  incident  cosmic  ray  particles  at 
the  LARC  location  (62°  12’  09”  S  -  58°  57’  42”  W,  South  Shetlands, 
King  George  Island,  Fildes  Bay  Ardley  Cove,  40  m  a.s.l.).  The  digital 
computer-based  studies  were  made  at  spaced  intervals  of  rigidity  in  a 
wide  range  extending  from  a  value  sufficiently  high  (20.00  GV),  such 
that  free  particle  access  is  ensured,  to  a  low  value  at  which  it  may 
confidently  be  assumed  that  access  would  be  forbidden  (2.00  GV)  for 
this  .Antarctic  location.  A  summary  of  the  obtained  results  is  given. 


1  Introduction 

Extraterrestrial  charged  particles  reaching  a  fixed  location  on  the  Earth 
travel  along  specific  paths  {trajectories)  through  the  geomagnetic  field.  The 
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knowledge  of  the  particle  direction  in  free  space  {asymptotic  direction),  i.e. 
prior  to  its  interaction  with  the  geomagnetic  field,  is  a  prerequisite  for  relat¬ 
ing  observed  cosmic-ray  intensity  variations  at  the  Earth  with  the  direction 
of  these  particles  in  the  near  Earth  environment  {e.g.  [1]  for  a  review).  In 
order  to  determine  which  energies  of  cosmic  ray  particles  are  able  to  reach  a 
selected  geographic  location,  it  is  necessary  to  perform  detailed  and  e.xten- 
sive  numerical  calculations  of  particle  trajectories  in  a  mathematical  model 
of  the  Earth  s  magnetic  field  (see  [2]  for  an  early  work).  In  fact,  the  Earth'.s 
magnetic  field  provides  a  shield  against  the  incoming  charged  particle  radi¬ 
ation.  the  particle  energy  threshold  is  high  in  the  Earth’s  equatorial  region, 
and  it  goes  towards  low  values  when  one  approaches  polar  latitudes. 

To  perform  trajectory  computations  independent  of  the  considered  parti¬ 
cle  species,  the  energy  (E)  is  replaced  by  the  rigidity  (R),  i.e.  the  momentum 
per  unit  charge.  The  name  cutoff  rigidity  is  used  to  indicate  the  minimum  R 
value  that  a  charged  particle  can  possess  to  arrive  from  a  specific  direction 
of  free  space  to  a  selected  measurement  point.  In  other  words,  the  cutoff 
rigidity  is  the  parameter  used  for  denoting  the  rigidity  value  below  which 
no  particles  can  be  detected  from  a  specific  direction  at  a  selected  location. 
This  value  is  dependent  on  zenith  and  azimuthal  directions  of  arrival.  Since 
most  of  cosmic  ray  particles  are  positively  charged,  there  exists  an  Ea$t-\Ve.sr 
effect  in  the  incoming  cosmic  ray  flux  to  a  detector,  with  higher  cutoff  rigid¬ 
ity  values  to  the  East,  and  lower  cutoff  rigidity  values  to  the  West.  Hence 
there  are  several  parameters  to  be  specified  before  exploiting  any  evaluation. 
They  are:  (i)  the  geographic  coordinates  of  the  detector  site;  (ii)  the  selected 
zenith  direcnon;  (iii)  the  chosen  azimuthal  angle;  (iv)  the  altitude  of  the  de¬ 
tector  location;  (v)  the  geomagnetic  status  (quiet  or  perturbed:  which  is 
often  described  by  pomagnetic  indices:  AE  for  high,  Kp  for  middle  and  Dst 
for  low  geomagnetic  latitudes)  at  the  measurement  time  and  the  baseline 
cjuiescent  geomagnetic  field  representation  for  the  appropriate  epoch. 

Calculations  should  be  initiated  at  a  R  value  well  above  the  maximum 
expected  East  cutoff,  and  continue,  at  discrete  intervals  S  =  AR  of  decrea.s- 
ing  R.  until  the  cutoff  is  certainly  surpassed.  A  careful  choice  of  the  R 
intervals  {i.e.  S)  is  needed,  because  for  many  experimental  .sites  (partic¬ 
ularly  those  at  middle  and  high  geomagnetic  coordinates)  the  transition 
of  allowed/forbidden  trajectories  is  not  unique.  The  name  cosmic  rny 
penumbra  is  used  for  denoting  the  region  between  the  first  (which  deter¬ 
mines  the  upper  rigidity  cutoff:  Ru)  and  the  last  (the  lower  rigidity  cutoff: 
Rl)  allowed/forbidden  transition  in  the  rigidity  space.  Moreover,  the  effec¬ 
tive  rigidity  cutoff  (Rc)  is  evaluated  by  subtracting  allowed  rigidities  in  the 
penumbra  from  the  Ru  value  (see  [3],  for  details  on  cutoff  rigidity  terminol- 
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^  ogy). 

\\  ithin  such  scenario  the  evaluation  of  induced  effects  on  cosmic  ray 
particles  by  the  long-term  changes  of  the  geomagnetic  field  at  the  Antarctic 
Laboratory  for  Cosmic  Rays  (LARC)  was  initiated.  The  first  results  were 
reported  by  Storini  et  al.  [4].  This  paper  presents  a  more  extensive  studv 
including  the  global  aspects  of  the  problem. 


2  Cosmic  ray  asymptotic  changes  for  LARC 

The  o-.vear  incremented  standard  models  for  the  quiescent  geomagnetic  field 
were  used  (from  Epoch  19o5.0/DGRF55  to  Epoch  1995.0/IGRG9.5)  for  eval¬ 
uating  cosmic  ray  induced  effects  by  the  secular  variation  of  the  geomagnetic 
field  at  the  LARC  (62^  12’  09”  S  -  58»  57'  42”  W,  South  Shetlands  Ki^ 
George  Island,  Fildes  Bay  Ardley  Cove,  40  m  a.s.l.)  location.  No  external 

T  “d  P)  for  an  attempt  in  using  the 

Tsvg89  -  Tsyganenko/1989  [6]  -  external  field). 

%  stated  our  trajectory  calculations  at  the  top  of  the  atmosphere  (as¬ 
sumed  20  km  high)  and.  working  backward,  to  traced  a  negative  proton  out 
through  the  geomagnetic  field  model.  This  corresponds  to  the  trajectorv 
0  a  positive  proton  coming  from  the  interplanetary  medium  through  the 
magnetosphere,  impacting  the  atmosphere  at  20  km  above  the  station  site 
and  then  creating  a  nuclear  cascade  down  to  the  detection  location  on  the 
xiaith  s  surfs^cs.  Th6  p3X3.m6t6rs  usGd  wore.’ 


(i)  geographic  coordinates:  62.20®  S  -  301.04®  E: 

(ii)  zenith  direction:  0®  {i.e.  vertical  direction); 

(iii)  azimuth  angle:  none; 

(iv)  altitude:  not  relevant; 


(v)  geomagnetic  status:  not  considered,  because  we  are  interested  on  effects 

induced  by  the  main  (internal)  field; 

(vi)  S  =  1.00  GV  for  20.00-10.00  GV,  S  =  0.10  GV  for  9.30-6.30  GV.  S  = 
0.0.5  GV  for  6.25-5.40  GV,  and  S  =  0.01  GV  for  5.39-2.00  GV. 

11  computations,  for  every  considered  epoch,  we  selected  the 

allowed  (particle  paths  going  to  the  interplanetary  space)  and  forbidden 
(particle  paths  intersecting  the  solid  Earth  or  paths  of  particles  which  have 


411 


and  1995  (righO^o^TThe' f«  1955  (left) 
.ndtcated  by  a  cross  surrounded  by  a  square  (as  in  the  oSeponeS  Igu^^'  '* 


“ ‘?,?s  /»•"  ‘  "T”  ■'*  ”  “C£ 

J  the  Earth  {i.e.  as  seen  fram  far  out  in  space)  for  the  two 

'‘t  rf.  P^n.l)%rL  .het  t 

c tear  that  until  the  considered  20  GV  particle  direction  (particle  rigidities 
increase  from  the  right  to  the  left  in  each  panel)  the  induced  change,  he 
aecular  variations  of  the  geomagnetic  field  are  significant  for  LARC  bcation 

Results  are  better  illustrated  in'  fig.  2,  where  the  particle  directions  in 

tennanced  toi  R  <  8  GV)  and  Southward  large  for  R  >  n  GV)  shifts 
are  evident  for  all  of  them.  Moreover,  the  longLdinal  spread  of  particte 
I  ’^JSidities  is  outstanding.  This  is  confirmed  at  4  GV  where 

e  particle  direction  shifted  from  the  Latin-American  sector  (19.55)  to  the 
.  Simi  one  (199o).  It  is  a  westward  shift  greater  than  240'  (see  fig  3) 

rhe^fpr.?  “>'”Ptotte  directions  of  the  incident  cosmic  ray  particles  at 
LARC  location  for  1995, 1990  and  1985  are  illustrated  in  fig.  4  for  rigidity 
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1955-1995  CHARGED-PARTICLE  ASYMPTOTIC  CHANGES  FOR  LARC 


Figure  2:  Particle  asymptotic  directions  at  LARC  location  from  1955  to  1995,  as 
Hg?dities'^^*”^  5-year  step  geomagnetic  field  models,  for  the  indicated  particle 


values  above  the  cosmic  ray  penumbra.  These  directions  are  again  projected 
on  a  map  of  the  Earth  (the  first  180°  of  longitude  are  repeated  on  the  right 
■  side  to  facilitate  particle  direction  identification).  Particle  directions  from  20 
GV  to  3.31  GV  (upper  panel),  3.42  GV  (middle  panel),  and  3.54  GV  (lower 
panel)  are  shown  as  circles  (see  the  left  and  right  sides  of  each  panel),  while 
a  circle  surrounded  by  a  square  denotes  the  directions  of  the  last  allowed 
GV  interval  (central  part  of  each  panel):  3.30  GV  -  3.21  GV  (upper  panel), 
3.41  GV  -  3.32  GV  (middle  panel)  and  3.53  GV  -  3.45  GV  (lower  panel)! 
The  Westward  shift  with  the  years  is  clearly  evident.  Changes  in  the  cosmic 
ray  penumbra  are  discussed  in  the  next  section. 

The  results  shown  above  for  vertically  incident  charged  particles  on 
L.ARC  detector  envisage  that  the  original  particle  direction  of  approach  (far 
away  from  the  Earth,  i.e.  the  asymptotic  direction)  can  be  very  different 
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1965^0^^5  fofThererSd^ection'MSLAMT*”" 

,he  variability  from“iVte;  »?he 


from  the  considered  detector  direction.  \everthelP9^  .  • 

“an- r 

E=~=SSS~ 

of  LARcic^T  “oasurements  (beginin. 

GLEfi^m  K  GLE51  (June  11  1991  m 

ofdbec  Iv  International  Da  a  (Demers 

c  ml  v“‘’  ''‘-"■®fi-"-ma3.i,;  (llj)  for  use  by  the  scie,  ifi 


3  LARC  geomagnetic  cutoff  or  threshold  rigidity 

Since  a  long  time  cosmic  ray  measurements  were  shown  to  be  affected  hv  tl„ 
secular  variation  of  the  geomagnetic  field.  Shea  and  SmLt  1101  aTj  •  • 

thirty- vear  interval  fl9‘?'5  tn  Aa  j  ^inahsing  a 

.  .  interval  (19J5  to  1965)  demonstrated  that  some  regions  of  the 
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IGRF95 


30*  60'  90’  126* ' i'so’  iso’  '210‘  '240*  276* '3o6*’330'  360’  ‘  30'' '  60*' '  90*'  ‘  1 20’ '  ISO' ' 


^DO-a-a, 


3.35  QV 


,.LDGRF90 


30'  60’  90*  120'  iso*  ISO’  210'  240*  arO’ '3o6’‘w  360’ '  30’' '  60*' '  90‘"i20'’ ISO* 


DGRF85 


30'  60'  90'  120'  150'  ISO'  '210*  240'  276' 'sOO'  ^'o'  SlBO' '  30'' '  60'' '  90'' '  120' '  IM' 

EASTERN  GEOGRAPHIC  LONGITUDE 

f  ’  directions  for  vertically  incident  cosmic  ray  particles  at  the 

LARC  location  as  calculated  using  the  DGRF85-DGRF90-IGRF95  geomagnetic 
field  models  (see  the  text  for  details). 


Earth,  particularly  the  Latin  American,  are  affected  by  significant  changes  in 
cosmic  ray  cutoff  rigidities  (see  also  [13,  14,  15],  among  others).  Cooper  and 
Simpson  [16].  using  the  Huancayo  neutron  monitor  data,  confirmed  such 
finding  for  Peru.  Figure  5  illustrates  the  changing  quiescent  geomagnetic 
field  from  1965  to  1995.  The  decrease  in  time  of  the  magnetic  field  intensity 
at  Huancayo  implies  an  increased  access  to  lower  rigidity  cosmic  rays  at  the 
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tX3RFI96S  •  GEOMAGNETIC  FIELD  INTENSITY  -  F  (nT) 


DGRFI975  -  GEOMAGNETIC  FIELD  INTENSITY  -  F  (nT) 


lor%t  F  at  the  surface  of  the  Ear, 


detection  point.  As  a  consequence,  an  increase  in  cosmic  ray  intensity  m 
the  years  is  expected  and  it  was  found.  Nowadays  the  Huancayo  detect 

P‘''“  cosmic  ray  observatories  axe  at  work  in  t 
Latin  American  sector  [17).  They  should  take  into  account  the  report' 

for  fhl  inlml'So  mo  components  at  LARC  locatic 

or  cne  interval  1955.0-2000.0,  using  one-year  increment.  They  were  derivf 

using  the  Synthesis  program  (version  3.0)  made  available  by  NOAA.  T1 
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dSit«^X®'Y™z?S‘?4Rr  f  ^  “?  ““PO”®""  (X.  Y,  Z)  and  their 

1QQ1  f j  (^ .  Y  ,  Z  )  at  LARC  location.  The  time  occurrence  of  1969,  1978  and 
1991  worldwide  geomagnetic  jerks  J^s,  j9i)  js  shown  in  lower  left  panel  (after 
[4];  arrows  indicate  the  corresponding  Y’  changes).  ' 


long  term  field  variability  at  LARC  is  clearly  emerging.  Figure  6  (lower  left 
panel)  shows  the  occurrence  of  fast  (<  5  years)  variabilities  {geomagnetic 
jerks)  in  the  slope  of  the  curve  of  the  annual  secular  variation  of  the  geo¬ 
magnetic  field  (e.^.  [18],  and  references  therein).  Their  effects  on  the  allowed 

cosmic  ray  rigidity  spectrum  at  a  measurement  point  should  be  further  in¬ 
vestigated. 

The  rigidity  spectrum  of  allowed  and  forbidden  particle  rigidities  in  the 
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ir(trpa„'el)  f^LARc'  As  It' 

penumbra  does  not  change  regularly  in  tlL  WWIe  R  d 

Steadily  from  1965  to  1995  the  R  vai  -n  ^  ^  decreases  nearly 

trend.  Certainly.  tL  fs  need  for  a  he»  “  T 

need  for  a  better  evaluation  of  possible  effects 
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Figure  8:  Results  from  Rc  evaluation  for  LARC  location  based  on  the  secular 
v^iation  of  the  quiescent  geomagnetic  field  (see  the  text  for  details).  The  data  fit 
of  the  1999  evaluation  [4]  is  shown  by  a  dashed  line. 


induced  by  the  South-American  geomagnetic  anomalies  and  by  the  world 
wide  geomagnetic  jerks  on  the  particle  paths  to  LARC.  Nevertheless,  we 
have  derived  from  them  the  effective  vertical  rigidity  cutoff  Rc  at  every 
epoch  and  they  are  reported  by  open  squares  in  fig.  8.  Recently,  Shea 
and  Smart  [19]  have  computed  vertical  effective  cutoff  rigidities,  from  1955 
to  1995,  for  a  long  list  of  neutron  monitor  stations.  LARC  location  was 
considered  again;  the  new  values  are  reported  by  filled  triangles  in  fig.  8. 
From  there,  it  is  possible  to  confirm  the  decreasing  trend  of  Rc  (dashed 
line)  at  LARC,  which  should  be  considered  in  any  long  term  study  of  the 
LARC  records. 

4  Conclusion 

Using  the  appropriate  Definitive  Geomagnetic  Reference  Field  coefficients 
for  epochs  1955  to  1990  (DGRF55  to  DGRF90)  and  the  provisional  ones  for 
1995  (IGRF95)  the  cosmic  ray  access  to  the  LARC  detector  was  studied  for 
quiescent  geomagnetic  conditions.  The  obtained  values  testify  a  decreasing 
cutoff  rigidity  trend  not  only  in  the  region  of  the  so-called  geomagnetic 
anomalies,  as  claimed  in  the  past  (see  Sect.  3),  but  also  Southward  of  it 
{i.e..  in  the  extreme  Latin  American  sector),  as  it  has  been  shown  for  the 
first  time  by  Storini  et  al.  [4].  The  reported  values  should  be  considered  as 
reference  values  for  additional  investigations. 
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